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ABSTRACT 


This  report  summarizes  a theoretical  investigation  of  the  current  problems 
of  nitramine  (composite)  propellant  combustion.  This  study  has,  as  its  distinc- 
tive feature,  a detailed  examination  of  the  condensed-phase  processes  in  the 
combustion  of  nitramine  propellants.  As  a consequence  of  a recently  developed 
model  for  the  combustion  of  ammonium  perchlorate  (AP)/composite  propellants, 
it  is  hypothesized  that  the  condensed-phase  degradation  of  the  nitramine  oxidizer 
particles  to  a vaporizable  state  is  the  overall  rate -limiting  step.  It  is  also 
assumed  that  the  gas-phase  details  are  secondary  in  importance  and  need  be 
studied  only  to  the  extent  of  supplying  the  correct  boundary  conditions  on  the 
condensed-phase/vapor-phase  heat  transfer.  Because  of  our  imprecise  under- 
standing of  the  gas-phase  processes  in  the  presence  of  combustion,  several 
plausible  models  are  considered  for  the  gas  phase.  It  is  found  that  all  of  the 
gas-phase  models  considered  lead  to  predictions  sufficiently  close  to  exper- 
imental trends  for  us  to  conclude  that  the  precise  details  of  gas-phase  processes 
are  not  of  critical  importance  in  determining  propellant  combustion  behavior. 
More  to  the  point,  we  are  led  to  believe  that  a thorough  examination  of  the  con- 
densed-phase details  may  be  sufficient  in  itself  not  only  to  interpret  most  of  the 
available  data  on  experimental  regression  rate  vs.  pressure  of  nitramine 
propellants,  but  also  to  aid  in  the  formulation  of  propellants  to  suit  our  needs. 
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SECTION  I 


INTRODUCTION 


The  limitations  that  the  highly  empirical  nature  of  solid  propellant  design 
place  on  its  technological  applications  are  revealed  clearly  whenever  new 
problems  are  encountered.  This  is  particularly  true  when  new  applications  are 
considered.  Recently,  the  nitramine /composite  propellants  have  been  recognized 
as  possessing  very  desirable  properties  for  certain  applications  (Ref.  1). 

Offsetting  the  thermodynamic  performance  superiority  is  the  combustion  behavior, 
which  has  shown  undesirable  slope  breaks  in  experimental  regression  rate  vs. 
pressure  curves,  The  undesirability  of  such  a characteristic  was  amply  described 
earlier  (Ref.  1).  It  was  thought  to  obviate  this  undesirable  behavior,  at  least  in 
the  operating  range  of  pressures  (approx,  below  30,  000  psi),  by  making  changes  in 
propellant  formulations.  This  was  to  be  done  in  such  a manner  as  not  to  adversely 
affect  the  good  features  of  performance,  etc.  In  the  absence  of  a theory  of  pro- 
pellant combustion,  however,  efforts  in  this  field  cannot  be  well  organized,  and 
: ti  the  absence  of  a working  model,  at  least,  such  efforts  are  conducted  almost 
totally  in  the  dark.  In  effect,  what  we  now  have  is  a set  of  experimental  data 
that  indicate  general  trends,  and  what  we  need,  as  a first  requirement,  is  a 
model  that  can  cohere  it l y interpret  most,  if  not  all,  of  the  available  data. 

A very  similar  situation  existed  in  the  related  field  of  AP/composite 
propellant  combustion  until  recently.  A host  of  superficially  diverse  data  could 
he  found  in  the  literature.  The  data  covered  the  degradation,  decomposition, 
and  sublimation  of  the  oxidizer  (AP),  and  the  degradation  of  the  binder  polymer; 
it  also  included  similar  data  on  AP  doped  with  various  chemicals,  hot  plate 
pyrolysis  data  on  AP  and  the  binder,  single  crystal  deflagration  data  on  AP, 
time-independent  data  on  burning  of  composites,  and  similar  time-independent 
data  on  composites  with  catalysts,  oscillatory  burning  of  composites,  and  all 
sorts  of  related  experiments.  An  attempt  was  made  to  interpret  this  mass  of 
data  coherently  in  terms  of  a few  simple  postulates  concerning  the  fundamental 
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processes  in  propellant  burning.  A model  was  presented  (Refs.  2-4)  under  the 
basic  hypothesis  that  the  fundamental  degradation  rate  of  AP  to  a vaporizable 
state  is  the  overall  rate-limiting  reaction  in  propellant  combustion.  This 
C1T/JPL  model  successfully  predicted  the  observed  trends,  including  a few 
that  were  inconsistent  with  previous  theoretical  predictions. 

In  the  present  study  this  model  has  been  explored  to  determine  its 
applicability  to  nitramine  propellant  combustion.  Aside  from  the  obvious 
difference  that  the  numerical  data  on  rate  constants  and  physicochemical  pro- 
perties would  be  different  from  those  for  AP/composites,  there  are  also  a 
few  subtle  differences,  and  these  are  mentioned  in  the  report.  In  general, 
however,  it  is  taken  as  a working  hypothesis  in  this  exploratory  study  that  the 
essential  details  of  the  model  are  applicable  to  nitramine  propellant  combustion. 
Justification  for  this  hypothesis  can  come  only  from  examining  the  results, 
although  - in  a weak  form  - an  indication  of  its  general  correctness  comes  from 
the  simple  reasoning  that  the  rate  - lim  iting  reactions  are  likely  to  be  in  the 
relatively  low  temperature,  condensed-phase  processes  for  nitramine  propellant 
combustion  also. 

Because  we  lack  some  of  the  input  data  needed  to  render  our  model 
completely  self-contained,  we  have  made  a parametric  study  of  the  problem. 

.'he  oxidizer  particle  size,  the  volumetric  loading,  and  the  melt  layer  thickness 
(explained  later)  are  varied  over  a range  in  model  propellant  formulations  to 
predict  theoretically  the  regression  rate  vs.  pressure  curves  and  related 
quantities  (surface  temperature,  gas  phase  combustion  zone  extent.  . .etc.). 

The  predicted  trends  are  found  to  be  close  to  the  experimentally  observed 
trends.  Most  of  the  numerical  constants  used  are  available  in  the  open  literature 
or  have  been  supplied  by  the  USAF.  The  remaining  few  parameters  are  not 
arbitrary.  They  have  been  the  subject  of  reasonable  estimates  in  the  past  and 
should  be  within  the*  grasp  of  experimental  efforts  in  the  near  future. 

The*  available  burning  rate  data  and  additional  US AF-supplied  data  are 
evaluated  in  Section  II,  after  a brief  literature  survey.  The  literature  survey  is 
restricted  to  those  models  of  propellant  burning  that  show  promise  of  application 
here.  The  effects  of  burning  rate  catalysts  and  oxidizer  particle  size  are  mentioned. 
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Suction  111  i*oi* s into  the  details  of  the  application  of  the  CIT/JPL  AP/ 
composite  propellant  model  to  nitramine  propellant  combustion.  This  material 
was  presented  in  the  monthly  technical  reports  made  by  JPL  to  the  U.S.  Air 
Face  over  the  period  of  this  contract.  The  assumptions  of  constant  wall 
temperature,  uniform  combustion  rate,  and  the  flame  standoff  distance  being 
equal  to  the  condensed  phase  heterogeneity  arc  used  to  develop  several  possible 
models  for  the  gas-phase  processes.  The  thickness  of  the  widely  accepted 
"melt  layer"  is  computed  using  a very  simple  uodel. 

In  general,  our  study  has  indicated  that  increasing  the  condensed-phase 
homogeneity  (and  hence  increasing  the  heat  transfer  rate  in  the  propellant)  or 
decreasing  the  oxidizer  particle  si2e  ought  to  have  beneficial  effects  in  reducing 
the  high  value  of  the  burning  rate  exponent  or  avoiding  the  slope -break  phenomena. 
These  can  be  regarded  as  design  hints  as  predicted  by  our  model.  Actual  quan- 
titative predictions  of  the  burning  rate  vs.  pressure  curves  for  nitramine/ 
composite  propellants  of  specified  formulations  will  require  information  on  a 
few  parameters  that  is  not  yet  available  and  can  or.'y  come  from  further  research. 
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SECTION  II 

EVALUATION  OF  AVAILABLE  LITERATURE  AND  DATA 


A.  AVAILABLE  LITERATURE 

A literature  survey  on  nitramine  combustion,  in  the  particular  context  of 
current  interest,  has  been  made  by  Ludwig  Stiefel  (Ref.  5).  Rather  than  repeat 
his  effort,  we  will  concentrate  on  those  papers  and  aspects  that  either  did  not 
receive  much  attention  in  Ref.  5 or  that  affect  our  own  work  directly.  The 
available  literature  can  logically  be  classified  in  two  categories  - as  experi- 
mental or  theoretical  work. 

I.  Experimental  Work 

Among  the  experimental  work  on  the  subject,  a large  number  have  to  do 
with  developmental  type  of  research  aimed  at  incr<  sing  our  understanding  of 
combustion  behavior  in  an  empirical  manner.  The  effects  of  small  variations  in 
catalysts,  in  binder  type  (energetic  or  inert),  and  in  particle  size,  have  all  been 
investigated  at  various  pressures.  The  available  data  should,  in  principle,  be 
helpful  in  designing  new  nitramine  propellants  for  specific  applications. 

However,  the  empirical  nature  of  all  of  these  efforts  (Refs.  6-9)  limits  their 
scope  where  direct  applicability  is  concerned.  Obviously  what  is  needed  is 
more  fundamental  experimental  work,  possibly  on  the  oxidizer  and  the  fuel 
separately,  even  if  such  efforts  appear  to  be  rather  isolated  from  technological 
applications.  It  is  encouraging  to  observe  that  two  such  efforts  have  indeed 
been  reported  in  the  open  literature. 

J. W.  Taylor  (Ref.  10)  found  in  his  work  on  PETN,  RDX,  and  HMX  that, 
for  any  given  pressure,  the  mass  burning  rate  was  practically  constant  at  all 
loading  densities  and  particle  sizes  of  pressed  strands  of  the  material  (up  to  the 
maximum  pressure  tested,  200  atm),  although  the  apparent  burning  rate  differed 
considerably  as  the  oxidizer  particle  size  was  varied  in  the  strands.  High-speed 
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motion  pictures  showed  unmistakable  evidence  of  a melt  layer  on  the  surface 
of  the  propellants.  He  postulated  that  the  burning  rate  is  controlled  by  the  melt 
layer  behavior  on  the  surface,  thus  accounting  for  the  negligible  influence  of 
particle  size  on  the  mass  burning  rate.  Taylor  also  conducted  experiments  on 
the  same  oxidizers,  PETN,  RDX,  and  HMX,  by  lightly  filling  tubes  of  paper 
and  Perspex  (PMMA)  in  various  mean  particle  sizes  (Ref.  11).  It  was  found 
that  beyond  a certain  pressure  the  linear  burning  rate  was  dramatically  higher, 
with  an  increase  in  the  pressure  exponent,  n (in  r = apn)  also.  He  attributed 
this  phenomenon  to  convective  heat  transfer  inside  the  pores.  This  "porous- 
bed"  burning,  although  very  relevant  to  the  overall  situation  in  the  cartridges, 
is  not  directly  relevant  to  burning  of  single  strands  of  nitramine  propellants, 
which  happens  to  be  our  prime  concern  in  this  report. 

The  work  of  Zimmer-Galler  (Ref.  1?.)  provides  us  with  the  numerical 
constants  for  the  degradation/decomposition  of  the  nitramines  of  interest  in 
high  pressure  applications.  Since  the  basic  rate  data  enter  into  all  fundamental 
models  in  a very  direct  manner,  the  work  of  Zimmer-Galler  has  proved 
extremely  useful  in  our  own  predictions.  An  interesting  feature  of  that  work 
(Ref.  12)  is  that  burning  rate  experiments  were  conducted  in  widely  different 
gas  environments  and  the  same  curve  of  regression  rate  vs.  pressure  was 
obtained.  This  indicates  that  the  fundamental  rate-limiting  reactions  are  unlikely 
to  be  in  the  gas  phase. 

An  examination  of  the  surfaces  of  the  quenched  samples  of  nitramine 
propellants  showed  a glazed  appearance,  leading  to  the  belief  that  a melt  layer 
existed  on  the  surface  during  combustion.  This  melt  layer,  w’hich  appears  to 
make  its  presence  felt  at  various  stages  of  nitramine  combustion  modeling,  is 
discussed  in  greater  detail  in  Section  III.  To  anticipate  some  ot  the  future 
discussion,  the  melt  layer  seems  to  figure  prominently  in  practically  all  of  the 
modeling  efforts,  but  whether  the  melt  layer  appears  before  or  after  the  break- 
point has  been  a matter  of  doubt.  At  any  rate,  the  kinetics  constants  determined 
by  Zimmer-Galler  have  been  consistently  used  throughout  our  work,  and  therein 
lies  the  value  of  that  work  for  our  discussions. 
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The  very  recent  experimental  work  on  RDX  decomposition  {Refs.  13  and 
14)  adds  nothing  fundamentally  new  to  the  field,  and  its  relevance  to  propellant 
combustion  studies  is  yet  to  be  demonstrated.  It  postulates  a mechanism  for 
the  decomposition  of  RDX  and,  as  such,  the  details  are  more  appropriately 
considered  theoretical  and  not  experimental. 

2.  Theoretical  Work 

The  amount  of  theoretical  work  on  nitramine  combustion  repoited  in  the 
literature  is  rather  limited  and  appears  to  have  been  done  only  in  the  general 
context  of  other  propellants.  Exceptions  do  exist  (Ref.  15),  but  are  few  in 
number.  Hence,  it  would  appear  that  only  those  theoretical  treatments  that 
are  sufficiently  general  to  include  in  their  scope  different  propellant  ingredients 
need  be  examined  in  some  detail.  The  basic  objective  of  such  an  examination 
would  be  to  extract  information  useful  in  arriving  at  ?.  rational  theory  of  nitra- 
mine combustion.  Unfortunately,  not  many  sufficiently  general  theories  on 
composite  propellant  combustion  exist.  The  well  known  theories  from  Princeton 
(Ref.  16)  admittedly  address  themselves  to  AP/composites  and  rely  rather 
heavily  on  the  available  AP  data.  Hence,  at  least  in  their  original  form,  the 
Princeton  theories  do  not  seem  to  have  direct  applicability  to  nitramine  combus- 
tion. 


The  modeling  efforts  at  Georgia  Institute  of  Technology  (Ref.  17)  are 
considering  sandwich  models  of  oxidizers  and  binders.  Such  studies  definitely 
add  to  our  knowledge  and  have  the  added  advantage  that  they  can  be  directly 
derived  by  simple  experiments,  some  of  which  have  already  been  performed  in 
AP/binder  systems  (Ref.  18).  It  is  not  clear,  however,  as  to  how  these  sand- 
which  studies  can  be  generalized  to  handle  the  considerably  more  complicated 
heterogeneous  mixture  of  oxidizer  and  binder  that  constitutes  a technology 
propellant  of  our  times.  Until  such  a method  of  applying  the  sandwich  results  to 
actual  propellants  is  indicated,  their  utility  with  regard  to  nitramine  combustion 
will  not  be  obvious. 
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That  essentially  leaves  us  with  the  BDP  (Beckstead..  Derr,  Price;  Refs. 
19-21)  model  of  composite  propellant  combustion,  which  is  gaining  popularity 
lately.  Here,  applicability  to  nitramine  combustion  is  definitely  possible; 
as  a matter  of  fact,  the  original  work  (Ref.  20)  actually  considered  such  a 
direct  application  to  HMX  (and  its  systems).  Extension  of  that  model  to  the 
specific  problem  on  hand  definitely  seems  feasible.  This  aspect  of  the  litera- 
ture will  not  be  considered  at  length  here  because  such  an  approach  is  the  main 
theme  of  a contemporary  study  at  the  original  BDP  center  (Lockheed  Propulsion 
Co.)  funded  by  AFOSR.  The  two  efforts  (LPC  and  JPL)  are  complementary, 
and  several  similarities  will  be  mentioned  below.  As  a matter  of  fact,  one  of 
the  main  objectives  of  the  present  work  is  to  explain  the  experimental  results 
of  the  LPC  work. 

The  CIT/JPL  model  of  composite  propellant  combustion  is  sufficiently 
general  to  be  applicable  to  a host  of  problems  in  composite  propellants.  In 
particular,  applications  to  current  problems  in  nitramine  propellant  combustion 
seem  feasible.  In  fact,  such  an  application  was  indeed  considered  and  the 
encouraging  results  obtained  were  presented  (Ref.  22)  at  the  Explosives  and 
Combustion  Meeting  held  at  the  Eglin  AFB  in  April  1974. 

At  present,  practically  all  of  the  useful  design  information  on  solid  rockets 
is  empirical.  Because  of  a general  lack  of  understanding  of  the  fundamentals, 
costly  programs  are  invariably  needed  in  any  solid  propellant  development. 

Hence  there  is  a need  for  work  at  a fundamental  level,  to  attempt  to  deempiricize 
propellant  combustion  data.  Within  the  framework  of  any  propellant  combustion 
nit  ode  1,  we  recognize  three  regimes  of  importance  to  propellant  combustion. 
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Only  the  essentials  of  the  CIT/JPL  model  are  presented  here;  additional  details 
are  contained  in  Refs.  2-4. 
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These  are  the  relatively  low-temperature  condensed  phase,  the  medium- 
temperature  vap  sr-phase /condensed-phase  interface  (wall),  and  the  relatively 
high-temperature  vapor  phase.  From  the  lact  that,  at  any  given  pressure,  we 
do  observe  a fairly  well-defined  regression  rate,  it  is  evident  that  an  overall 
rate-limiting  reaction  is  in  operation  in  the  system.  The  basic  idea  behind  the 
CIT/JPL  model  is  that  the  overall  rate-limiting  reaction  in  propellant  combus- 
tion is  likely  to  be  in  the  relatively  low-temperature  condensed  phase.  The 
vapor-phase  heat  transfer  to  the  propellant  surface  matches  the  requirements 
of  the  degradation  reactions  and  the  sensible  enthalpy  rise.  Since  many  of  the 
oxidizers  used  in  propellants  are  similar  - crystalline  solid  particles  and  the 
binders  are  polymers  of  one  form  or  another  - and  since  the  basic  processes 
in  the  condensed  phase  (degradation  and  decomposition)  are  also  similar,  it 
would  seem  worthwhile  to  investigate  the  possibility  of  a single  model  being 
applicable  to  all  such  propellant  combustior. . If  the  key  processes  behind  the 
overall  burning  rate  are  in  the  vapor  phase,  the  prospects  of  the  success  of 
such  an  investigation  appear  bleak  because  r,f  the  rather  varied  compositions 
that  occur  in  the  vapor  phase  as  determined  by  propellant  chemistry. 

Since  the  CIT/JPL  model  was  based  on  the  fundamental  assumption  that 
the  key  processes  are  in  the  condensed  phase,  fairly  detailed  calculations  were 
performed  of  the  regression  rate  vs.  pressure  curves  and  the  related  quantities 
(flame  standoff  distance  from  the  surface,  wall  temperature  variations,  etc.). 
The  theoretical  results  concerned  themselves  only  with  AP  and  AP/composite 
propellants.  The  predictions  were  found  to  be  close  to  experimental  observa- 
tions in  various  applications  including  hot  plats-  pyrolysis  data  of  various  inves- 
tigators (compiled  by  Powling,  Ref.  23),  AP  single  crystal  deflagration  data, 
initial  temperature  sensitivity,  and  the  time  independent  combustion  data  of 
propellants.  In  addition,  the  same  model  also  ave  predictions  on  oscillatory 
combustion  characteristics  (response  functions)  that  were  close  to  experimental 
observations,  but  unexplained  by  previous  theories  (Ref.  24).  Those  gas -phase 
details  needed  to  supply  the  proper  boundary  conditions  on  the  condensed  phase 
are  modeled  to  be  consistent  with  two  of  the  popu.ar  pictures.  It  is  found  that 
either  is  capable  of  handling  the  situation  adequately,  thereby  de -emphasizing 
the  importance  of  the  precise  details  of  the  gas-phase  processes  to  propellant 
c ombustion. 
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B.  USAF-SUPPLIED  BURNING  RATE  DATA 


The  USAF- supplied  burning  rate  data  on  nitramines  was  evaluated  in  the 
framework  of  the  CIT/JPL  model.  Rather  a large  number  of  parameters  were 
varied,  and  the  experimental  output  also  varied;  the  parameters  of  direct 
relevance  and  importance  to  the  CiT/JPL  model  are  as  follows:  the  burning 
rate  exponents  before  and  after  the  break  point,  the  propellant  formulation, 
the  particle  size  of  the  oxidizer  and,  as  explained  later,  the  ratio  of  the  thermal 
diffusivity  to  oxidizer  particle  size.  This  ratio  (K  / a)  has  the  dimensions  of 
regression  rate  and  figures  prominently  in  our  theory  of  nitramine  combustion. 
The  USAF-supplied  data  on  regression  rate  vs.  pressure  is  extensive  and 
could  be  discussed  at  length  in  the  context  of  our  work.  However,  it  is  con- 
sidered more  appropriate  to  prepare  a tabLe  of  the  relevant  quantities  to  enable 
the  reader  to  ascertain  the  general  trends  in  the  data  (Table  1).  As  can  be 
seen,  although  the  pressure  at  the  breakpoint  is  considerably  different  for 
different  propellant  formulations,  the  regression  rates  at  the  breakpoints  do  not 
differ  greatly.  In  general,  this  leads  us  to  believe  that  the  basic  hypotheses 
behind  our  condensed  phase  model  may  be  valid,  .since,  using  only  condensed- 
phase  variables  (thermal  diffusivity,  particle  size,  and  regression  rate),  we 
have  been  able  to  unify  some  of  the  superficially  diverse  data  on  nitramine 
combustion. 

C.  CATALYSTS  AND  THEIR  EFFECTS 

There  does  not  appear  to  exist  in  the  open  literature  a systematic  study  of 
the  nitramine  propellant  behavior  as  catalysts  are  varied  in  type  or  concentration. 
As  already  discussed  the  available  information  is  almost  entirely  empirical  and 
does  not  lend  itself  to  a systematic  classification  of  the  effects  of  catalysts  on 
burning  rate  or  other  effects.  This  state  of  affairs  with  catalysts  is  recognized 
by  Stiefel  also  (Ref.  5). 

Specifically,  there  does  not  appear  to  be  any  evidence  (Ref.  I)  that  the 
addition  of  catalysts  would  either  eliminate  or  make  a significant  improvement 
in  this  burning-rate  slope-break  phenomenon. 
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D.  OXIDIZER  PARTICLE  SIZE  EFFECTS 


As  already  mentioned,  the  effect  of  oxidizer  particle  size  offers  a promising 
area  for  further  investigation,  since  variations  in  this  physical  parameter  appear 
to  affect  the  slope-break  phenomenon  rather  strongly  (Ref.  1).  In  general, 
decreasing  the  oxidizer  particle  size  seems  to  help  avoid  both  the  high  slope  and 
the  slope-break  phenomenon,  as  is  shown  in  Fig.  1.  In  fact,  using  extremely 
small  oxidizer  particles  ( 1 p)  seems  to  remove  the  slope-break  phenomenon 
completely  from  the  operating  range  of  pressures  (Ref.  15).  The  experimental 
data  available  in  this  respect  are  summarized  in  Table  1.  As  can  be  seen, 
there  is  no  conclusive  evidence  as  to  the  precise  nature  of  the  systematic 
variations  in  the  slope-break  variations  with  oxidizer  particle  size.  However, 
there  apparently  exist  clear  indications  of  these  beneficial  effects  in  other  ways 
(Refs.  1 and  15).  A systematic  investigation  of  the  effects  of  decreasing  oxidizer 
particle  size,  intended  to  provide  a solid  ground  for  working  models  on  this  known 
beneficial  effect,  is  in  progress  at  the  Eglin  AFB  and  the  experimental  data  are 
being  received  periodically  at  JPL.  In  the  near  future,  it  should  be  possible  to 
process  the  systematic  data  in  terms  of  the  CIT/JPL  model. 
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Figure  1.  Effect  of  RDX  Particle  Sice  on  Measured  Regression  Kate 
Versus  Pressure,  Replotted  From  Reference  12 
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SECTION  III 


POSTULATED  MECHANISMS  OF  HIGH  PRESSURE  DECOMPOSITION 
A.  BASIC  POSTULATES 

The  combustion  of  composite  propellants  involves  extremely  complicated 
processes  in  all  of  the  three  regions  of  interest  (condensed  phase,  wall  inter- 
face, and  vapor  phase).  Any  completely  realistic  model  would  be  so  interactive 
analytically  that  it  would  be  impractical  to  handle.  What  is  needed  is  a model 
for  the  key  processes  that  retains  a degree  of  similarity  to  the  actual  physics 
and  chemistry,  at  least  to  the  extent  of  representing  the  fundamental  parameters 
and  phenomena.  Since  an  adequate  theoretical  understanding  of  most  of  the 
details  of  propellant  burning  does  not  exist,  the  only  real  test  of  a model  is  the 
agreement  (or  otherwise)  of  the  predicted  results  with  experimentally  observed 
trends.  However,  good  agreement  in  one  regime  is  no  indication  of  the  merits 
of  the  model  when  extrapolations  to  other  regions  are  considered.  If  a model 
consistently  predicts  results  close  to  experimental  observations  in  several 
different  regimes,  it  is  tempting  to  apply  the  model  to  other  systems  as  well. 

The  basic  concept  behind  the  C1T/JPL  model  is  shown  schematically  in 
Fig.  2.  The  fundamental  rate -limiting  reaction  is  hypothesized  to  be  the 
degradation  of  the  oxidizer  crystals  in  a thin  melt  layer  on  the  interface  between 
the  oxidizer  and  the  binder.  Hence,  the  geometry  of  the  propellant,  as  deter- 
mined by  the  formulation,  can  be  related  to  the  degradation/decomposition  of 
the  propellant.  It  is  to  be  clearly  remembered,  as  has  been  emphasized 
several  times  (Refs.  3 and  4),  that  the  model  is  an  average  representation 
only  and  is  NOT  to  be  interpreted  literally. 

Concentrating  our  attention  on  the  condensed  phase  only,  we  may  write 
the  one -dimensional  energy  equation  applied  between  the  deep  solid  (oo)  and  the 
vapor/solid  interface  wall  (o). 

k + per  4^  - D pc  B exp  (-E/RT) 
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Figure  2.  Propellant  Model 


The  source  (or  sink)  term  on  the  right-hand  side  of  the  energy  equation  contains 
terms  related  to  the  actual  degree  of  degradation  at  each  plane  in  the  propellant. 
This  can  be  related  to  the  fundamental  degradation  equation,  which  is,  in  the 
particular  case  at  hand,  an  Arrhenius  form  of  the  equation.  The  pro -exponential 
constant  B (in  k = B exp  (-E/RT))  is  assumed  to  be  linearly  dependent  on  the 
chamber  pressure.  The  justification  for  this  pressure  dependence  comes  from 
the  reasoning  that  the  actual  degradation  is  brought  about  by  the  diffusion  of  a 
catalytic  species  in  the  melt  layer  saown.  The  equations  applied  to  AP / 
composite  propellants  lead  to  the-  regression  rate  expression: 
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Detailed  information  concerning  this  equation  is  given  in  Refs.  2-4.  To  briefly 
recapitulate  the  definitions  of  the  symbols, 

K is  the  thermal  diffusivity  of  the  propellant. 

'I'  is  the  thickness  of  the  interfacial  melt  layer  on  the  oxidizer  particles. 

v is  the  volume  fraction  of  the  oxidizer  particles  in  the  propellant. 

is  the  reference  value  of  the  pre -exponential  constant  in  the  degradatio 

h is  defined  as  D/c  (Tw-T  ),  where  D is  the  heat  required  by  unit  mass 
of  the  oxidizer  to  be  degraded  to  the  fundamental  individual  repeating 
units. 

c is  the  specific  heat  of  the  propellant  (solid), 

Tq  is  the  temperature  of  the  unaffected  propellant. 

FSV  is  the  parameter  Fragment  Size  Vaporizing,  which  essentially 

quantifies  the  extent  of  propellant  degradation  at  the  vaporization 
step. 

As  can  be  seen,  "P,  the  interfacial  melt  layer  thickness,  is  the  only 
parameter  not  precisely  determined  at  the  present  time.  However,  since 
reasonable  estimates  of  its  values  can  be  made,  it  is  not  a free  parameter. 

For  example,  a that  is  5-10%  of  the  oxidizer  particle  size  seems  reasonable 
while  40-50%  certainly  does  not.  The  numerical  value  of  the  wall  temperature 
T , has  to  be  determined  through  proper  matching  with  the  gas-phase  details, 
and  several  different  models  have  been  tried  for  the  gas  phase.  Actually, 
experimentally  measured  values  of  the  wall  temperature  may  also  be  used 
when  available.  All  of  these  approaches  were  used  to  predict  various  quantitie 
of  interest  in  the  combustion  of  AP/composite  propellants  (Refs.  2-41. 

Because  of  these  and  other  (Ref.  26)  successes  of  the  preliminary  work, 
the  CIT/JPL  mode!  was  applied  in  the  present  study  to  the  combustion  of 
nitramino  propellants.  The  basic  aim  was  to  postulate  mechanisms  of 
degradation/decomposition,  so  as  to  remove  from  the  operating  range  the 
discrete  slope  breaks  in  the  experimental  regression  rate  vs.  pressure  curves 
Essentially,  the  same  analytical  picture  as  for  AP  propellants  was  used. 
Obviously,  the  numerical  constants  were  appropriately  different.  Also,  in 
the  original  CIT/JPL  model,  as  applied  to  the  combustion  of  AP/composite 
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propellants,  the  concept  of  Fragment  Size  Vaporizing  (F;7V)  was  introduced  to 
signify  the  extent  of  propellant  degradation  at  the  vaporization  step.  An  effort 
was  made  to  determine  the  numerical  value  of  FSV  independently  of  propellant 
combustion  details  by  observing  similarities  of  vapor  pressure  equilibrium 
data  of  hydrocarbons.  Here,  in  the  combustion  of  HMX,  RDX,  etc.  , this 
second  aspect  of  FSV  relation  seems  to  need  some  modifications.  First,  the 
pressures  of  interest  are  so  high  that,  at  all  reasonable  values  of  wall  tempera- 
tures, the  FSV  is  far  less  than  unity,  which  makes  degradation  of  the  funda- 
mental unit  in  the  oxidizer  RDX  and  HMX  seem  to  be  important.  The  FSV  rule, 
in  its  original  form  (Ref.  2),  is  not  applicable  to  such  a case.  Second,  the 
fundamental  unit  in  the  crystal  (HMX  and  RDX)  is  not  small  or  simple  (as  in  AP) 
and  therefore  the  concept  of  FSV  needs  direful  interpretation. 


Because  of  the  considerations  discussed  above,  the  experimentally 
measured  values  of  the  degradation  rates  were  used  directly.  These  values 
are,  inherently,  influenced  by  variations  of  the  surface  details  during  combustion. 
Basically,  vvc  are  assuming  that  the  experimentally  determined  fundamental 
degradation  rate  constants  are  directly  applicable  at  pressures  higher  than 
those  at  which  they  were  measured.  The  validity  of  such  an  approach  can  be 
established  only  by  examining  the  results.  The  results  obtained  in  the  present 
study  do  appear  reasonable  in  quantitative  comparison  with  experimental  data. 
Thus,  it  would  appear  that  the  new  assumption  concerning  the  vaporization 
step,  in  generalizing  the  original  C.1T/JPL  model  to  nitramines,  is  a valid 
approach.  Besides,  as  can  be  seen,  the  uncertainties  in  this  respect  (i.e.  , 
those  concerning  the  value  of  FSV)  are  likely  to  make  only  minor  variations  in 
the  quantitative  regression  rate  predictions  and  hardly  anv  at  all  in  the  quali- 
tative trends.  Thus,  even  if  the  new  assumption  on  1\S\  should  prove  wrong 
later,  our  results  here  would  still  retain  their  utility. 


As  will  become  clear  later,  the  basic  as*-  option  behind  the  derived 
equation  (Eq.  (2)|  is  that  the  condensed  phase  n at  rial  (the  propellant)  may  be 
considered  homogeneous  for  the  purpose  of  heat  transfer  calculations.  It  is 
this  assumption  that  enables  us  to  write  the  energy  equation  in  the  first  place. 
When  this  assumption  breaks  down,  the  result  |Eq.  (2)1  cannot  lie  valid. 
Obviously,  a second  model  is  needed  to  represent  the  propellant  combustion 
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mechanics.  The  new  model  is  shown  in  Fig.  3.  It  is  assumed  that  all  of  the 
degradation  reactions  now  take  place  in  a surface  layer  of  thickness,  i. 


Turning  our  attention  to  the  elusive  gas-phase  details,  the  popular  flame 
sheet  model  is  envisioned  as  a possible  means  of  exploring  the  wall  temperature 
variations  [and  regression  rates  (Ref.  2),  too].  The  essence  of  the  flame  sheet 
approximation  is  that  all  of  the  gas-phase  reactions  are  confined  to  a thin  zone 
parallel  to  the  surface  but  displaced  a distance  (X*)  from  it.  In  other  cases,  it 
is  possible  that  the  reactions  in  the  gas  phase  are  uniformly  distributed  and  are 
characterized  by  a mass  consumption  rate  of  m"'  = A ■ P,  where  A is  a constant 
and  P the  pressure.  With  this  model,  we  may  predict  the  wall  .emperature  and 
regression  rate  variations  if  the  value  of  A is  known. 
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Figure  3.  High  Pressure  Behavior 
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The  task  ahead  of  us  now  becomes  fairly  clear.  The  trends  in  propellant 
combustion  behavior  are  to  be  examined  as  some  of  these  parameters  are 
varied.  Consideration  is  also  to  be  given  to  the  applicability  to  nitramine 
propellant  combustion  and  the  implications  of  the  results  obtained  in  light  of  the 
available  experimental  data.  Such  an  approach  was  undertaken,  and  the 
periodic  progress  was  summarized  in  monthly  reports.  The  results  are 
presented  here. 

B.  EFFECT  OF  VARIATIONS  IN  WALL  TEMPERATURE  ON  THE 

SLOPE  OF  THE  BURNING  RATE  LAW  CURVE 

Breaks  in  the  slope  of  solid  propellant  burning  rate  vs.  pressure  relation- 
ships are  generally  indicative  of  discrete  dimensional  effects.  Discrete 
dimensions  in  heterogeneous  propellants  are  obviously  associated  with  the 
oxidizer  particle  [it  is  recalled  that  no  slope  breaks  have  been  reported  in  tests 
with  liquid  oxidizers  such  as  iso-DMED  (Ref.  7)].  Naturally,  the  question  is 
"under  what  conditions  do  the  discrete  dimensions  manifest  themselves?"  The 
natural  dimension  in  propellant  burning  is  the  propellant  thermal  depth,  K /T 
(propellant  t'  <"”-mal  diffusivity  t linear  regression  rate).  Wrhen  this  dimension 
is  large  compared  to  the  oxidizer  particle  size,  a,  we  expect  discrete  dimen- 
sional affects  to  be  smeared  out.  On  the  contrary,  when  k/~  is  small  compared 
to  a,  discrete  dimensional  effects  should  be  evident. 

In  the  present  section,  the  breakpoints  are  predicted  under  the  assumption 
that  they  occur  when  the  a racteriotic  thermal  depth  in  the  condensed  phase 
becomes  comparahle  to  or  less  than  the  characteristic  heterogeneity  scale  in 
the  propellant  (oxidizer  particle  size).  For  tne  pro -breakpoint  region,  the 
effects  of  oxidizer  particle  size  on  the  linear  regression  : ate  were  not  included 
at  this  stage  of  the  study.  These  effects  are  available  in  the  term  (b  v/a), 
which  is  assumed  here  to  be  constant,  i.e.  , at  a specified  oxidizer  volumetric 
loading,  v,  the  ratio,  ^/a,  of  the  interfacial  melt  layer  thickness,  to  the 
oxidizer  particle  size,  a,  is  assumed  to  be  constant  at  a reasonable  value  of 
5%.  Explicit  variations  of  the  oxidizer  particle  size  throughout  the  entire 
regime  are  considered  later.  At  the  present  stage  of  this  discussion,  the 
explicit  variations  in  the  post- breakpoint  regime  are  considered.  At  the  point 


of  slope  discontinuity,  i.  e.  , the  transition  from  subsurface  reaction  rate 
control  to  surface  reaction  rate  control,  the  matching  of  the  two  regression 
rates  determines  the  value  of  the  only  external  parameter  (the  melt  layer 
thickness,  t ) which  is  not  arbitrary. 

C.  CALCULATIONS  ASSUMING  CONSTANT  WALL  TEMPERATURE 

In  the  calculations  of  the  present  section  constant  wall  temperatures 
have  been  assumed.  Of  course,  the  correct  way  of  treating  the  problem 
is  to  make  no  assumption  regarding  the  wall  temperature,  but  rather  to  allow 
it  (Tw)  to  be  self-determined  by  the  gas-phase  details;  that  is  the  next  step, 
reserved  for  a later  section.  However,  we  may  make  a few  comments  here  on 
the  possible  implications  of  such  a treatment.  It  is  known  that  the  use  of  the 
gas -phase  reaction  rate  expression  in  its  simplest  form,  m,u  = A • P,  results 
in  a pressure  exponent  slightly  higher  than  that  given  by  the  constant  Tw  model. 
That  is  to  say,  the  values  of  the  exponent,  n,  revealed  by  the  present  calcu- 
lations would  increase  slightly  both  before  and  after  the  breakpoint.  To  hazard 
a guess,  the  n's  may  increase  from  the  present  values  of  0.  5 and  1 to  a 0.  6 
and  « 1.2,  which  are  much  closer  to  typic  1 experimental  values  (Thiokol  and 
Hercules).  The  main  point  to  note  is  that  the  present  section  gives  a specific 
calculation  only.  More  sophisticated  calculations  wall  be  described  in  sub- 
sequent sections. 

Slight  modifications  have  been  made  in  the  basic  CIT/.TPL  theory  (Refs.  2 
and  3)  in  attempting  to  apply  it  to  the  vastly  different  nitrarnine  gun  propellant 
system.  The  concept  of  FSV  does  not  appear  to  be  directly  applicable: 

1)  The  pressures  are  so  high  that,  at  all  reasonable  values  of  wall 
temperatures,  the  FSV  is  far  less  than  unity. 

2)  Since  the  fundamental  unit  in  the  crystal  is  not  small  or  simple  for 
either  RDX  or  HMX,  the  FSV  concept  needs  careful  interpretation. 

In  light  of  the  above,  the  experimentally  measured  values  of  the  degra- 
dation rates  will  be  used  directly.  These  values  are,  inherently,  influenced  by 


variations  of  the  surface  details  in  degradation.  The  details  of  the  procedure 
employed  are*. 


In  place  of  gas-phase  reaction  rates  (to  be  used  in  subsection  III-D), 
reasonable  values  of  the  wall  temperatures  are  assumed  (600°K 
and  650°K). 

A modified  form  of  subsurface  degradation  equation  is  used  to 
predict  the  linear  regression  rate. 

Recognizing  that  the  subsurface -reaction  model  is  of  limited 
validity  when  the  condensed  phase  is  not  homogeneous  (i.  e,  , when 
k! r < a,  the  predictions  are  stopped  at  r = k! a.) 

At  the  point  of  discontinuity  (in  slope),  the  surface  reaction  rate 
equation  is  used  and  the  prediction  of  r vs.  P is  continued.  There 
are  no  arbitrary  parameters  at  this  stage;  the  value  of  the  only 
external  parameter,  the  melt  layer  thickness,  i , is  determined  by 
matching  the  regression  rates  at  the  breakpoint 


The  numerical  values  used  are  given  in  Table  2. 


TABLE  2.  NUMERICAL  VALUES  USED  (RDX)  (CONSTANT  WALL  TEMPERATURE) 


Parameter 


The  equations  used  are  as  follows: 


1)  ‘Before  the  breakpoint: 


,,  -E/RT 

A BPe 
a o 


(E/RTw) 


T - T 
w o 


[P  in  atm  ] 


2)  At  the  breakpoint  r' : 

r“  = which  yields  values  of  “r  as  follows: 


a in  p 

2 

10 

30 

50 

100 

A , / 

r m cm/sec 

7.  5 

1.  5 

0.  5 

0.  3 

0.  15 

r1  in  in.  /-ec 

2.  95 

0.  59 

0.  196 

0.  118 

0.  059 
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REGRESSION 


D. 


CALCULATIONS  ASSUMING  UNI  FORM  COMBUSTION  IN 
THE  GAS  PHASE 


In  this  section  we  extend  the  previous  work  to  include  consideration  of 
gas-phase  reactions,  instead  of  assuming  a constant  wall  temperature.  In 
the  process,  we  describe  the  preparation  of  a computer  program;  later  in  this 
report  improvements  to  this  program  will  be  described.  The  program  is 
sufficiently  general  to  permit  determination  of  the  regression  rate  as  a function 
of  pressure,  when  the  condensed-phase  and  gas -phase  parameters  are  specified. 

1.  Method 

In  accordance  with  the  analyses  in  Ref.  2,  the  regression  rate  is 
determined  by  the  matching  of  the  gas -phase  mass  combustion  rate  to  the 
condensed-phase  mass  generation  rate.  The  value  of  the  wall  temperature 
enters  the  expressions  for  both  of  these  rates  and  iterative  procedures  are 
needed  for  solution,  since  transcendental  functions  are  involved.  The  general 
method  is  described  in  detail  on  pp.  35-38  (particularly  on  page  37)  of  Ref.  2. 

In  our  application  the  method  is  modified  by  elimination  of  the  parameter  FSV, 
which  is  believed  to  be  of  questionable  applicability  to  the  high  pressure  com- 
bustion of  complex  substances  like  nitramines. 

Hence,  the  following  equations  (which  result  from  a simplification  of  the 
general  equations  prese  ited  in  Ref.  2)  are  used: 

1)  Condensed-phase  mass  generation  rate: 
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2)  Gas-phase  mass  consumption  rate  (nondimcn sional ): 


A = 


K . Q ■ m1 1 1 

p . c . (T  - T ) r' 
r p w o 


(4) 


where  m'"  is  the  uniform  combustion  rate  in  the  gas  phase  (gm  cm 
sec  S and  Q is  the  heat  of  combustion  (cal  gm  ^). 


Matching  the  heat  flux  from  the  gas  to  the  condensed  phase  leads  to 


(A  • 0 + C3  exp  (C)  + C4  - 1 = T i - 1 4 h (5) 

' L \v  L o 

where  C = In  Z 

z = - a/c3 

C3  = I + h - A 
C4  = 1 - C3 

h = nondimensional  heat  of  degradation,  D/c  (T  - T ) 

It  is  evident  that  there  are  no  free  parameters  in  the  system.  The 
gas -phase  reaction  rate,  fn’",  uniquely  specifies  the  regression  rate  since  the 
other  variables  in  the  system  {Bq,  K»  D,  Q,  c,  pi  Tq1  all  have  very  definite 
numerical  values. 

The  gas -phase  mass  consumption  is  now  taken  as  adequately  represented 
by  two  parameters  A and  m 


tii"'  = A Pm  (6) 

This  is  an  assumption  in  the  theory.  However,  it  has  been  widely  used 
and  is  seen,  even  in  the  simplest  case  (m  5 1),  to  yield  very  reasonable 
results  for  AP /composite  propellants.  The  main  support  for  the  above  form 
of  gas -phase  reaction  rate,  aside  from  the  fact  that  it  has  been  highly  successful 
in  other  application,  is  described  in  detail  in  Rei.  2 (pp.  ^2-^1. 
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Since  the  mass  consumption  rate  in  the  gas  phase  is  assumed  to  be 
completely  determined  by  inte rmolecula t diffusion  of  the  reactants  in  the 
propellant,  the  numerical  values  of  the  constants  are  unlikely  to  be  influenced 
by  the  chemical  nature  of  those  vapors. 

Such  reasoning  leads,  as  a first  approximation,  to  the  same  numerical 
values  of  A and  m as  were  used  in  other  propellant  applications.  This  question 
of  the  numerical  values  of  A and  m has  to  be  eventually  settled  through  actual 
experimental  measurements  of  the  gas-phase  reaction  rate.  The  situation  is 
very  analogous  to  the  condensed-phase  parameters  E and  Bq,  which  come  from 
more  fundamental  experiments  (DTA.  TCA,  DSC,  etc.  ). 

2.  Computer  Program 


NOTE 

A listing  of  this  computer  program  ano 
the  deck  of  cards  were  sent  to  Eglin  AFB 
earlier.  The  listing  is  self-explanatory. 

Briefly,  the  computer  program  has  in  it  the  values  of  Bq,  E,  k,  c,  (j, 

and  T . It  reads  in  values  of  A,  m,  D,  Q,  trial  Tw,  ATW»  a,  and  v.  The 

above  values  form  one  set  in  a data  card.  The  number  of  such  sets  fthe  number 

of  data  cards)  is  specified  in  the  program  by  the  variable  NUMBER.  The 

program  then  uses  the  trial  value  of  Tw  to  compute  independently  the  two 

sides  of  the  equation  (Eq.  5)  (called  W1  and  W2  in  the  program).  The  two 

independently  calculated  values  must  be  equal  for  a proper  solution.  A 5% 

error  bound  has  been  established  in  the  program.  A 5%  error  in  Eq.  (5) 

results  usually  in  a 1 or  2%  error  bound  on  the  regression  rate,  which  is  felt 

to  be  adequate  at  this  stage.  The  error  bound  can  be  easily  altered.  If  the 

two  values  W1  and  W2  do  not  match  at  the  initial  trial  value  of  T , one  of  two 

w 

things  can  happen: 

1)  If  the  wall  temperature  is  too  LOW,  it  is  gradually  incremented 

in  steps  of  ATw  (called  DT  in  the  program)  until  a favorable  agree- 
ment is  obtained.  In  the  process,  if  it  is  found  necessary,  the 
incremental  step  size  ATw  (DT),  is  made  smaller  by  factors  of  10 
at  a time. 


2)  If  the  initial  wall  temperature  is  too  HIGH,  the  following  message 
should  appear: 

THE  INITIAL  VALUE  OF  TRIAL  WALL  TEMPERATURE, 
TRIAL  T,  IS  TOO  HIGH.  CHOOSE  A LOWER  VALUE. 

The  regression  rate  is  computed  for  four  values  of  pressure:  1,000, 

11,000,  21,000,  and  31,000  psi.  The  output  also  contains  the  values  of  A and 
m used,  the  wall  temperature  computed,  and  the  value  of  the  regression  rate 
at  which  breaks  are  likely  to  occur  in  the  standard  regression  rate  vs.  pressure 
plot.  For  one  set  of  the  parameters  the  results  obtained  are  plotted  in  Fig.  6. 

The  details  of  ^/a,  v,  binder,  etc.  arc  not  pursued  at  this  stage.  Th  s 
point  was  discussed  in  the  previous  section. 

It  is  noted  that,  at  this  stage,  the  program  is  preliminary  in  form.  It 
has  been  found  to  work  for  a few  sets  of  values,  for  all  of  which  m - 1.  It  has 
not  been  tested  extensively  for  various  values  of  the  parameters  or  to  optimize 
the  cost  considerations. 
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Figure  6.  Calculated  Regression  Rate  Versus  Pressure  Assuming 
Uniform  Combustion  in  the  Gas  Phase  (m"*  * A-P) 

2f» 


E. 


PARAMETRIC  VARIATIONS  OF  ii>  AND  a 


In  this  section  the  previous  work  is  extended  to  include  cases  in  which  the 
thickness  of  the  interfacial  melt  layer  is  not  constant,  but  variable,  and  in  the 
process  a computer  program  is  prepared  that  is  an  improvement  over  the 
previous  effort  (subsection  III-D-2).  In  addition,  the  other  important  parameter, 
namely  the  oxidizer  particle  size,  is  varied  over  the  range  of  general  interest. 
The  program  is  simple  enough  to  easily  handle  variations  in  the  volumetric 
loading  of  the  oxidizer  in  the  propellant. 

In  the  basic  equation,  Eq.  (2),  all  of  the  following  arc  determined  by  the 
propellant  formulation:  the  volumetric  loading  of  the  oxidizer  in  the  propellant, 
v;  the  oxidizer  particle  size,  a;  the  thermal  diffusivity,  K\  and  the  kinetics 
constants  Bq  and  E.  The  initial  temperature,  T , and  the  chamber  pressure,  P, 
are  determined  by  the  experimental  conditions;  the  propellant  wall  temperature, 
Tw>  is  determined  by  the  gas-phase  energetics  and  fluid  dynamics;  and  the 
intcrfacial  melt  layer  thickness,  has  to  come  from  experimental  measure- 
ments. Tw  and  & are  the  only  two  variables  that  do  not  have  unambiguous 
values  at  this  stage  of  the  analysis.  In  Section  III-D-2,  the  value  of  Tw  was 
determined  for  two  different  values  of  the  gas-phase  reaction  rate  parameters. 

It  was  found  that  the  Tw  did  not  vary  with  pressure  when  the  gas-phase  reaction 
rate  was  assumed  to  be  controlled  by  molecular  mixing  of  the  fuel  and  oxidizer 
vapors.  Based  on  the  above  exploratory  study  of  the  possible  variations  in  Tw> 
it  was  concluded  that  the  assumption  of  constant  T is  reasonably  justified  for 
nitramine  propellants  of  given  formulations.  It  is  not  evident  that  the  wall 
temperature  i^  constant,  only  that  this  may  be  a reasonably  good  assumption 
at  this  stage.  Of  course,  this  point  has  to  be  settled  through  careful  measure- 
ments of  the  gas-phase  reaction  rates.  Unfortunately,  such  experimental  data 
arc  unavailable  at  this  time;  they  are  being  obtained  in  a separate  project  at 
JPL,  In  any  case,  the  wall  temperature  is  assumed  not  to  vary  with  pressure 
in  the  present  calculations,  although  it  is  recognized  that  the  wall  temperature 
may  vary  with  variations  in  other  formulation  parameters  (such  as  the  oxidizer 
particle  size  for  example).  In  fact,  the  results  obtained  in  the  present  section 
strongly  suggest  that  the  wall  temperature  probably  decreases  with  increasing 
oxidizer  particle  size. 
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Using  the  above  background,  the  only  parameter  whose  value  is  not 
thoroughly  settled  at  this  stage,  the  melt  layer  thickness,  ip,  is  varied 
parametrically.  The  basic  idea  is  to  have  available  results  that  will  find  ready 
application  once  measurements  are  made  of  ip  variations.  It  has  been  previously 
assumed  that  ip  was  5%  of  a.  Now  its  value  is  varied  from  1%  to  9%  to  see  its 
effect  on  the  breakpoint  in  the  regression  rate  vs.  pressure  curves.  In  Fq.  (2) 
it  is  seen  that  the  regression  rate  varies  as  the  square  root  of  ft.  Even 
though  the  variation  is  mild,  it  can  nevertheless  have  a noticeable  effect 
on  the  breakpoint.  A subtle  point  is  that,  although  increases  in  oxidizer 
particle  size  decrease  the  regression  rate,  this  result  is  valid  only  when  all 
the  other  parameters  are  held  fixed.  For  example,  if  the  interfacial  melt  layer 
thickness,  ip,  is  increased  at  the  same  time,  the  regression  rate  can  actually 
increase  with  increases  in  particle  size.  Consequently,  when  we  specify  ip/a 
as  a certain  percentage,  the  results  have  to  be  examined  with  care. 

The  breakpoint  regression  rate  is  predicted,  as  above,  by 

r = if/a  (7) 

At  regression  rate  values  exceeding  *r,  since  the  propellant  cannot  be  considered 
homogeneous  anymore,  surface  reactions  in  a molt  layer  of  thickness  f are 
postulated  and  the  regression  rate  is  now  given  by 


r = f R P exp  (-F/RT  ) (8) 

o w 

The  value  of  the  melt  layer  thickness  at  the  wall  is  not  arbitrary,  but  is 
determined  by  proper  matching  a_t  the  breakpoint. 

The  procedure  is  as  follows: 

1)  The  breakpoint,  “r,  is  determined  by  r = */a 
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j*^,'****|**l^''*‘''r ii^i*' . 


T 


The  breakpoint  pressure  is  determined  hy 


r2  (E/RT  ) 

w T 

£ = _ 

*'  V exp  <-E/RTw> 


The  melt  layer  thickness  at  the  breakpoint  is  determined  by 


B • P . exp  (-E/RT  I 
o ^ w 


It  is  then  assumed  that  the  melt  layer  thickness  is  constant  in  the 
post-breakpoint  regime,  at  least  in  a reasonable  range  of  pressure 
variations. 

Computer  Program 


NOTE 

A listing  and  printout  of  this  computer 
program  were  sent  to  Eglin  AFB  earlier. 

The  listing  should  be  self-explanatory. 

THERDF  is  the  thermal  diffusivity,  K 
EL  is  the  surface  melt  layer  thickness,  I 
PS1  is  the  interfacial  melt  layer  thickness,  C 


The  computer  program  reads  in  values  of  the  wall  temperature,  T , and 
volumetric  loading  of  the  oxidizer,  v,  in  the  propellant.  It  computes  the  linear 
regression  at  five  discrete  values  of  the  pressure:  1,000,  11,000,  21,000, 

11,000,  and  41,000  psia.  The  program  also  considers  four  values  of  the  oxidizer 
particle  size  (2  p,  10  p,  50  p,  and  100  p)  and  five  values  of  the  interfacial  melt 
layer  thickness  (l?j,  3To,  5%,  7%,  and  9%  of  the  oxidizer  particle  diameter). 

The  program  automatically  selects  the  correct  regression  rate  equation, 
depending  on  the  location,  as  pre -breakpoint  or  post-breakpoint.  This  is  an 
improvement  over  the  previous  efforts.  In  the  pre -breakpoint  regime  the 


r 
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surface  melt  layer  19  nonexistent  (at  least  within  the  framework  of  Refs.  2 
and  3),  and  hence  the  melt  layer  thickness,  is  printed  out  as  0.0.  The  rest  of 
the  results  should  be  self-evident.  The  computer  costs  are  quite  low  ($1.06). 

2.  Discussion 

The  results  show  that  increases  in  the  melt  layer  thickness  increase 
the  regression  rate,  as  anticipated.  In  Refs.  2 and  3 the  formula  [Eq.  (2)] 
was  derived  under  the  assumption  that  p is  small  compared  to  a.  Hence,  when 
ti)/a  approaches  or  exceeds  C,  1 (for  example),  the  results  may  not  be  valid.  In 
such  a case,  physically  speaking,  the  interfacial  melt  layer  thickness  on  the 
oxidizer  particle  is  becoming  comparable  to  the  oxidizer  particle  diameter 
itself  and  the  basic  model  becomes  of  questionable  applicability. 

In  comparison  with  experimental  values  of  linear  regression  rates,  the 
present  results  indicate  that  the  wall  temperature,  T^,  probably  decreases 
with  increases  in  oxidizer  particle  size.  As  can  be  seen  in  the  printout,  the 
linear  regression  rate  predictions  are  reasonable  at  the  higher  values  of  the 
wall  tempi  rature  (~  630°K)  and  the  higher  values  of  the  molt  layer  thickness 
u 9%)  for  the  smalle r -oxidize .* -pa r tic le  propellant,  while  the  regression  rates 
of  the  larger -oxidizer -particle  propellants  seem  reasonable  at  the  lower  values 
of  the  wall  temperature.  Such  a variation  in  the  wall  temperature  is  reasonable 
within  the  framework  of  a wide  variety  of  gas-phase  models  (G1)F  models,  flame 
sheet  models,  uniform  combustion  rate  models,  etc.).  This  is  so  because,  as 
oxidizer  particle  size  increases,  the  gas -phase  combustion  zone  moves  farther 
and  farther  away  from  the  propellant  surface  (wall).  This  results  in  a decrease 
in  wall  temperature,  since  the  basic  energetics  are  not  affected.  However, 
this  picture  needs  careful  interpretation  in  the  post-breakpoint  regime  since,  hy 
the  very  nature  of  the  surface  reaction  model,  oxidizer  particle  size  is  no  longer 
a significant  factor. 

F.  GAS -PHASE  DETAILS  AS  INFLUENCED  BY  T11E  OXIDl/.KK 
PARTICLE  SIZE 

Although  we  have  all  along  been  using  a condensed -phase  theorv  that 
hypothesizes  that  most  of  the  interesting  details  in  propellant  combustion  are 
controlled  by  the  condensed -phase  processes,  we  da  need  to  consider  a few 
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essential  details  of  the  gas-phase  processes  in  order  to  obtain  a self-contained 
solution.  In  their  simplest  form  the  gas-phase  details  are  adequately  contained 
in  the  value  of  the  temperature  at  the  gas /solid  interface  popularly  called  the 
"wall  temperature,  T . " The  value  of  the  wall  temperature  determines  the 
linear  regression  rate  of  the  propellant  through  the  condensed-phase  equation 
lEq.  (2 >].  In  our  previous  sections  two  seemingly  different  models  were 
proposed  for  the  gas-phase  processes;  both  were  based  on  the  assumption  that 
the  molecular  mixing  rate  of  fuel  and  oxidizer  vapors  controls  the  rate  of 
chemical  reaction  in  the  gas  phase: 

1)  The  assumption  of  constant  wall  temperature,  with  a flame  sheet 
model  for  the  gas -phase  combustion. 

2)  The  assumption  of  uniform  combustion  in  the  gas  phase. 

Neither  of  the  above  two  models  attempted  to  interpret  the  gas -phase 
details  in  a mechanistic  way.  The  present  section  mal<es  a first  attempt  to 
overcome  some  of  the  limitations  and  actually  relate  the  gas -phase  details  to 
propellant  formulations. 
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1. 


Method 
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The  basic  ideas  of  the  analytical  approach  of  the  present  section  are 
described  below.  The  flame  zone,  or  the  zone  of  vigorous  combustion,  is 
established  over  the  propellant  surface.  Because  of  the  discrete  nature  of  gas 
evolution  from  the  particles,  it  is  expected  that  considerable  inhomogen'  ity 
exists  in  the  gas  phase  above  the  burning  surface.  In  fact,  such  gas -phase 
inhomogenities  are  clearly  visible  in  the  photographs  of  AP/composite  propellant 
combustion  (Ref.  27).  Far  downstream,  however,  such  inhomogenities  dis- 
appear because  the  combustion  reactions  reach  completion.  It  is  clear,  there- 
fore, that  the  scale  of  gas-phase  inhomogeneity  (reflected  in  the  flame  zone 
dimensions)  is  related  to  the  physical  extent  of  condensed-phase  heterogeneity 
(reflected  in  the  oxidizer  particle  size).  Hence,  a quantitative  relation  is 
sought  between  these  two  variables.  The  idealized  representation  of  the  flame 
zone  dimensions  is  the  "flame  standoff  distance."  Analytically,  on-  can  predict 
the  propellant  regression  rate  if  the  flame  standoff  distance  is  specified.  The 
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flame  standoff  distance  not  only  varies  with  the  condensed-phase  heterogeneity, 
but  also  with  the  chamber  pressure.  Any  theoretical  modeling  of  the  flame 
standoff  distance  must  predict  both  of  these  variations.  At  the  present  stage 
we  concentrate  only  on  the  variations  with  the  oxidizer  particle  size.  The 
essence  of  the  present  approach  is  written  as 

X*  oc  a (11) 

or  X;':  = (Constant)  • a 

A linear  dependence  was  temporarily  assumed  between  the  flame  standoff 
distance  and  the  oxidizer  particle  size.  Based  on  physical  considerations,  it 
can  be  expected  that  larger  oxidizer  pa ' tides  would  result  in  the  flame  zone  being 
established  farther  away  from  the  surface.  However,  the  assumed  1 1 ne a r 
dependence  [Eq.  (11)]  is  by  no  means  obvious.  It  is  merely  one  of  several 
possible  assumptions  that  enable  us  to  obtain  analytical  solutions  at  this  stage. 
This  form  was  chosen  purely  for  the  simplicity  it  affords.  As  can  be  readily 
appreciated,  the  representation  of  Eq.  (11)  needs  modification  when  different 
pressures  are  considered.  For  example,  the  numerical  value  of  the  constant 
in  Eq.  (11)  can  be  expressed  as  a function  of  pressure  to  incorporate  the  fact 
that,  for  a given  propellant  formulation  (oxidizer  particle  size),  the  flame 
standoff  distance  decreases  with  increases  in  pressure. 

The  analytical  solution  is  straightforward.  The  propellant  formulation 
specifies  the  values  of  a and  v,  along  with  the  standard  values  of  Bq,  E, 

K (sk/pc),  and  T . A reasonable  (constant)  value  of  the  interfacial  melt  layer 
thickness,  l li,  is  assumed  (=  5%  of  a).  The  value  of  the  wall  temperature 
uniquely  determines  the  regression  rate.  However,  the  wall  temperature 
cannot  have  an  ambiguous  value  if  it  is  required  that  the  adiabatic  flame 
temperature  (determined  by  the  propellant  formulation)  be  reached  exactly  at 
the  "flame  standoff  distance,"  as  specified  by  Eq.  (11).  Hence,  the  gas -phase 
energy  equation  applied  between  the  wall*  and  the  flame  sheet  is  solved  for 
the  heat  transfer  rate  at  the  wall  (wall  temperature  gradient).  There  exists 
only  one  value  of  the  wall  temperature  at  whii  h the  heat  transfer  rate  into  the 
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condensed  phase  from  the  gas  phase  vaporizes  the  propellant  material  at 
exactly  the  rate  determined  by  the  condensed -phase  chemical  kinetic  degra- 
dation equation  [Eq.  (2)]. 

2,  Computer  Program 


NOTE 

The  listing  of  this  computer  program  was 
sent  to  Eglin  AFB  earlier.  The  listing 
should  be  self-explanatory,  in  view  of  the 
programs  sent  before. 

The  computer  program  is  written  such  that  the  value  of  the  temperature 
at  the  flame  sheet  is  printed  out  for  incremental  values  of  the  wall  temperature, 
Tw#  A knowledge  of  the  adiabatic  flame  temperature,  T^,  for  the  propellant 
thus  enables  us  to  locate  the  exact  value  of  the  linear  regression  rate  from  the 
computer  printout.  The  wall  temperature  is  also  printed  out  along  with  the 
regression  rate.  This  provides  an  independent  "check"  on  the  solution.  Since 
the  wall  temperature  range  for  the  nitramine  propellant  combustion  is  fairly 
well  known,  the  wall  temperature  actually  obtained  (for  the  correct  flame 
temperature)  can  be  useful  in  evaluating  the  procedure. 

3.  Numerical  Data 


The  RDX  constants  in  Table  3 were  used  for  the  representative  propellant 
considered  here. 

4.  Results 

The  results  are  presented  in  Figs.  7 and  8.  It  is  most  encouraging  that: 

1)  The  wall  temperatures  are  predicted  to  be  about  600eK,  which 
is  the  generally  accepted  value  for  nitramine  propellant 
combustion. 
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TABLE  3.  NUMERICAL  VALUES  USED 


Parameter 

Val  ue 

Source 

Activation  energy  of  degradation,  E 

48,  000  cal/mole 

Ref.  12 

Pre -exponential  constant,  B 

irJ9.1  -1 

10  sec 

Ref.  12 

Thermal  diffusivity,  K = k/pc 

-4  2 

15x10  cm  /sec 

Eglin  AFB  data 

Volumetric  loading,  v,  of  RDX  in 
the  inert  binder  propellant 

80% 

Typical  value 

Ratio  of  interfacial  melt  layer 
thickness  to  oxidizer  particle 
size,  </)/a 

5% 

Assumed  as 
reasonable 

Oxidizer  particle  size 

1 0 p - 1 00  p in 
steps  of  10  p 

Flame  temperature  of  mono- 
propellant RDX 

3282  °K 

Ref.  12 

Flame  temperature  (adiabatic) 
of  propellant  with  20%  inert 
binde  r 

2640°K 

Eglin  AFB  data 

Constant  in  X*  = (constant)  • a 

1 

Assumed 

WTTtt 


WALL  TEMPERATURE,  ‘K 

>40  600  660 


0.?  0.4  0.6  0.8  1.0  1.2 


LINEAR  REGRESSION  RATE,  cm/sec 


Figure  7.  Calculated  Regression  Rate,  Propellant  Wall  Temperature,  and 
Oxidizer  Particle  Size  Relationships  Assuming  A Simple  Flame  Standoff 
Distance  Gas-Phase  Model;  Pressure  = 1000  psia 
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OXIDIZER  PARTICLE  SIZE, 

Figure  8.  Calculated  Regression  Rate,  Propellant  Wall  Temperature,  and 
Oxidizer  Particle  Size  Relationships  Assuming  A Simple  Flame  Standoff 
Distance  Gas-Phase  Model;  Pressure  =*  21,000  psia 


2) 


The  regression  rates  are  predicted  to  be  around  0.4  cm/sec  at 
1000  psia,  which  is  again  reasonable. 

3)  The  wall  temperature  and  the  regression  rates  are  predicted  to 
decrease  with  increases  in  the  oxidizer  particle  size. 

The  last  of  these  is  the  mo3t  significant  result  of  the  present  program. 
The  widely  observed  experimental  trend  of  decreasing  regression  rate  with 
increases  in  the  oxidizer  particle  size  is  correctly  predicted  by  this  extremely 
simple  modeling.  This  particular  version  of  the  model  does  not  have  for  its 
aim  the  prediction  of  X*  variations  with  pressure.  However,  an  analysis  in 
subsection  III-H  which  predicts  the  X*  variations  with  mean  pressure  will  show 
that  the  X*  decreases  nearly  logarithmically  with  pressure  increase.  When 
that  fact  is  incorporated  into  the  present  model,  the  simple  model  should  be 
remarkably  useful  indeed. 

5.  Breakpoints 

The  breakpoints  are  predicted  within  the  general  context  of  the  CIT/JPL 
model  as  occurrmg  when  the  characteristic  thermal  depth  in  the  solid  becomes 
less  than  the  characteristic  heterogeneity  scale  (oxidizer  particle  size)  in  the 
propellant.  As  reported  in  the  September  1 973  technical  monthly  report  of  this 
program,  the  breakpoints  are  predicted  to  be  as  follows: 
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G.  MELT  LAYER  COMPUTATIONS 


It  is  suggested  that  the  observed  melting  phenomenon  is  associated  with 
the  thermal  gradient  in  the  condensed  phase.  The  wall  temperature  is 
generally  believed  to  be  approximately  constant,  at  about  600  - 650”K.  The 
melting  points  of  RDX  and  KMX  are  reported  (Eglin  AF13  data)  to  be  204°C 
and  284°C,  respectively.  Hence,  the  propellant  material,  in  its  travel  from  the 
low’  temperature  environment  ( ~ 300°K)  to  the  high  temperature  wall,  may 
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undergo  melting  after  the  melting  point;  that  is,  the  temperature  difference 

(Tw  - is  traversed  through  a melt  layer.  It  is  recognized  that  the  argument 

is  not  rigorously  valid,  since  the  quoted  melting  points  refer  to  an  equilibrium 

phenomenon,  whereas  propellant  combustion  is  a nonequilibrium  phenomenon, 

_ 3 

with  residence  times  of  the  order  of  10  sec.  However,  this  crude  argument 
is  thought  to  be  adequate  as  a first  description  of  a complex  process.  Neglecting, 
for  the  moment,  the  heat  of  degradation  in  the  condensed  phase,  the  solution  to 
the  energy  equation. 


is 


+ 


0 


-rx 

K 


with  x measured  from  the  wall  into  the  propellant  material.  Thus,  at  any 
regression  rate  (and  hence  pressure),  the  melting  point,  specifies  tl 

thickness,  £,  of  the  melt  layer  for  a specified  wall  temperature,  T 


w 


K / T - T \ 

I = — In  ( -.j,  W “ rnelt  layer  thickness  (12) 

\ mp  " o) 

In  the  very  short  time  available  (typically  in  the  region  T — f/r  ~10  ^ sec)  it  is 
unlikely  that  mixing  of  binder  and  oxidizer  can  take  place  thoroughly  through 
the  melt  layer;  consequently,  the  present  model,  which  considers  spherical 
oxidizer  particles  with  a thin  interfacial  melt  layer  on  them,  may  still  be  valid. 
As  has  been  pointed  out  in  earlier  reports  (Ref.  5,  for  example),  the  CIT/JPL 
model  is  intended  to  be  a highly  idealized  (but  useful)  representation  of  the 
complex  mechanics  in  propellant  combustion. 


The  data  included  in  the  September  1973  monthly  technical  report  to 
Eglin  AFB  is  the  basis  for  the  present  computations.  As  a typical  example,  the 
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r vs.  P curve  for  1 0 (j,  oxidizer  particles  and  600°K  wall  temperature  are 
chosen.  The  following  table  is  easily  generated  through 


T - T -rl 
= e— 

rn  rr^ 

1 w o 


(204  + 271)  - 100 
600  - 300 


15x10' 


microns 


Pressure,  psia 


1300 

2000 

3000 

3650 

5000 

6000 

7000 

8000 

| 10,000 

0.  3 

0.  37 

0.  45 

0.  5 

0.  58 

0.  64 

0.  69 

0.  74 

0.  825 

26.  4 

21.4 

17.  6 

15.  8 

l 3.  6 

12.  4 

11.  5 

10.  7 
l 

9.6 

The  results  are  plotted  in  Fig.  9. 

The  following  points  should  be  borne  in  mind: 

1)  Heat  of  degradation  has  not  been  considered.  The  wall 

temperature  gradients  in  actuality  will  be  rteeper  and,  consequently, 
the  melt  layers  will  be  thinner  than  predicted,  as  shown  in  Fig.  10. 


Wall  temperatures  are  not  precisely  known  at  this  time  and  600  °K 
is  an  approximate  number  only.  A higher  value  of  T^  leads  to  a 
lower  value  of  the  melt  laver  thickness. 

Binder  interactions  have  been  totally  ignored.  These  will  lower 
the  melt  layer  thickness  because  the  binder  heat  sink  leads  to  a 
steeper  wall  temperature. 


4) 


Variations  in  the  value  of  the  thermal  diffusivity  with  either  Tw 
or  a could  alter  the  values  presented  in  Fig.  9. 

5)  The  computations  were  performed  using  thermal  and  chemical 

parameters  (T^  , B,  c,  p,  . , . ) for  a typical  RDX  propellant. 
HMX  has  a higher  melting  point- -284  °C.  Hence,  the  melt  layers 
are  likely  to  be  substantially  thinner  than  for  RDX  propellants. 


It  is  suggested  that  the  melt  layer  thickness  computed  here  be  regarded 
as  an  upper  bound  only.  As  shown  in  Fig.  9,  the  theoretical  results  cannot 
be  extended  to  the  possible  "disappearance"  of  the  melt  layer  at  higher  pressures 
because  of  limitations  on  the  validity  of  Eq.  (12)  used  to  compute  the  thickness. 
Equation  (12)  assumes  that  the  propellant  material  may  be  considered  homo- 
geneous for  the  purposes  of  heat  transfer  calculations.  The  homogeneous  solid 
assumption  is  of  questionable  validity  beyond  the  regression  rate  at  which  the 
characteristic  thermal  depth,  tf/r,  equals  the  characteristic  heterogeneity 
scale,  a. 


H.  VARIATIONS  OF  FLAME  STANDOFF  DISTANCE  WITH  PRESSURE 


Here  we  consider  the  flame  standoff  distance  variations  in  greater  detail. 
The  flame  standoff  distance  is  an  analytically  convenient  concept  which  enables 
us  to  write  a simple  expression  for  the  heat  transfer  rate  from  the  vapor  phase 
to  the  condensed  phase.  Obviously,  for  a given  energetic  scheme  (as  determined 
by  the  propellant  formulation),  this  heat  transfer  rate  must  depend  on  the 
chamber  pressures.  This  is  nothing  more  than  a statement  of  the  fact  that  the 
regression  rate  is  pressure  dependent. 


A theoretical  prediction  of  the  variations  in  the  flame  standoff  distance 
with  pressure  would  enable  us  to  further  evaluate  our  model  regarding  its 
applicability  to  nitraminc  propellant  combustion.  Also,  as  explained  later, 
quantitative  experimental  determination  of  the  flame  standoff  distance  with 
pressure,  which  is  relatively  easy  to  measure  (compared  with  the  wall  tempera- 
ture, Tw,  or  the  gas -phase  combustion  rate,  m"’),  may  actually  help  us 
determine  some  of  the  other  unknowns. 
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The  computations  here  are  subject  to  the  same  limitations  as  before. 

For  example,  the  breakpoints  are  still  predicted  as  being  given  by  the  ratio 
of  the  thermal  diffusivity  to  the  oxidizer  particle  size.  Also,  predictions  based 
on  the  condensed-phase  model  are  not  vax’d  beyond  the  breakpoint. 

The  key  elements  of  the  computation  in  this  section  are  as  follows; 

1)  The  variations  in  the  flame  standoff  distance  are  c nnputod  based 
on  a reference  value  at  a reference  state.  It  is  assumed  that  the 
flame  standoff  distance  is  proportional  to  the  condensed-phase 
heterogeneity  (the  oxidizer  particle  size).  Moreover,  at  the 
reference  pressure  of  1000  psia,  the  flame  standoff  distance  is 
equated  to  twice  the  oxidizer  particle  size.  This  equality  is 
somewhat  arbitrar-y , but  is  not  thought  to  be  of  critical  importance 
to  the  general  conclusions  of  this  study. 

2)  The  wall  temperature  is  assumed  to  remain  constant  as  the  pressure 
is  varied.  However,  Tw  variations  with  oxidizer  particle  size  are 
considered. 

The  computational  procedure  is  straightforward.  The  propellant  formu- 
lation specifies: 

1)  The  oxidizer  particle  size,  a. 

2)  The  volumetric  loading  of  the  oxidizer,  v. 

3)  The  flame  temperature,  T^>  which  is  assumed  not  to  vary 
appreciably  with  pressure. 

We  assume  a reasonable  value  of  the  interfacial  melt  layer  thickness 
(5%  of  the  oxidizer  particle  size,  a).  The  linear  regression  rate  of  the 
propellant  is  computed  at  incremental  values  of  the  wall  temperature,  Tw» 
using  Eq.  (2),  The  flame  standoff  distance  is  computed  through  X*  = C^9^r» 
and  the  melt  layer  thickness  (Fig.  11)  is  computed  using  the  procedure  dis- 
cussed in  the  previous  section  (Eq.  (12)].  The  computations  have  been  carried 
out  for  wide  variations  in  the  wall  temperature  for  oxidizer  particle  sizes  of 
20,  30,  40,  50,  and  60  p and  at  pressures  of  1000,  11,000,  21,000,  31,000, 
and  41 , 000  psia. 
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Figure  11.  Calculated  Variation  in  Surface  Melt  Layer  Thickness  Kith 
Pressure  and  Oxidizer  Particle  Sire 


There  exist  several  different  ways  of  interpreting  the  extensive  data 
obtained.  Figure  12  presents  the  variations  in  the  flame  standoff  distance  with 
pressure  variations.  The  calculations  realistically  predict  that  the  flame 
standoff  distance,  X - , decreases  with  pressure.  It  must  be  borne  in  mind  that 
these  predictions,  using  Fq.  (2)  are  not  valid  beyond  the  breakpoint,  r.  How- 
ever, in  view  of  the  crudeness  of  the  breakpoint  predictions  (orde r-of-magnitude 
validity  only),  the  predicted  line  are  continued  with  broken  lines  up  to  the 
highest  pressure,  41,000  nsia.  The  variations  in  the  oxidizer  particle  size, 
as  they  influence  the  flame  standoff  distance,  are  considered.  It  is  clear  that 
the  predictions  in  Fig.  12  are  probably  accentuating  this  effect  somewhat.  It 
is  known  from  experiment  that  a variation  in  the  oxidizer  particle  size  from  ^0 
to  20  p.  is  unlikely  to  increase  the  regression  rate  by  more  than  a factor  of  2 
However,  the  trends  are  predicted  well. 
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figure  12.  Calculated  Variation  in  Fla»e  Standoff  Distance  and 
Regression  Rate  With  Pressure 
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SECTION  IV 


EXAMINATION  OF  THE  MODEL  IN  THE  CONTEXT  OF  THE 
AFOSR/AFATL  EXPLOSIVE  COMBUSTION  MEETING 


In  this  section  the  authors  would  Like  to  examine  the  CIT/JPL  theoretical 
model  for  nitramine  propellant  combustion  in  light  of  the  information  gained 
from  the  meeting  presentations  and  discussions.  Specifically,  it  is  intended 
to  show  that  the  CIT/JPL  model,  though  superficially  restrictive,  is  sufficiently 
general  to  be  consistent  with  current  ideas  on  nitramine  combustion. 

A.  KEY  ELEMENTS 

1)  Gas-Phase  Processes:  In  general,  many  investigators  seem  to  be 

amenable  to  considering  the  condensed  phase  as  the  site  of  the  key 
processes.  The  arguments  against  the  importance  of  the  gas-phase 
details  have  for  their  supporting  evidence  the  experiments  of 
/.imme  r-Galler  (Ref.  Id),  where  essentially  the  same  regression 
rate  behavior  was  obtained  in  widely  different  gas  environments. 

The  motion  pictures  of  Cohen  (Ref.  15)  seem  to  show  that  there 
exists  a fairly  clearly  defined  bright  "flame"  zone  some  distance 
( ~ 1 00  p)  from  the  propellant  surface.  This  may  be  associated 
with  the  flame  sheet  model  that  was  used  in  subsection  III-H. 

d)  Melt  Layer:  There  appears  to  be  almost  universal  agreement  that, 

at  least  in  the  pre-breakpoint  region,  a melt  Layer  exists  on  the 
surface  of  nitramine  propellants.  Motion  pictures  jVisnov  (Ref.  28), 
for  example],  surface  examination  of  quenched  samples  (Zimmer- 
Galler  (Ref.  Id),  for  example],  and  scanning  electron  microscope 
photographs  (Ref.  29),  all  support  the  concept  of  a melt  layer. 

There  also  appears  to  be  an  almost  equally  prevalent  belief  that 
this  melt  layer  gets  thinner  as  the  pressure  increases. 


3)  Post-Breakpoint  Regime:  It.  was  the  general  belief  that  this  regime 
of  combustion  is  not  well  understood.  More  experimental  data  are 
awaited.  Even  analytical  models  tend  to  be  speculative  without 
reliable  experimental  data. 

B.  DISCUSSION 

In  analogy  with  the  available  studies  on  AP  single  crystal  deflagration 
[see  for  example,  Guirao  and  Williams  (Ref.  30)],  it  is  tempting  to  postulate 
the  existence  of  a melt  layer  that  gets  progressively  thinner  with  increasing 
chamber  pressure.  (Actually,  the  regression  rate  is  more  relevant  than  the 
pressure.)  Beyond  a certain  value  of  She  regression  rale,  the  melt  layer 
disappears  completely,  and  hence  the  regression  ran-  vs.  pressure  curve  is 
qualitatively  different.  Such  ideas  on  AP  deflagration  need  careful  interpretation 
before  they  can  be  applied  to  nitramine  combustion.  It  is  to  be  remembered 
that  the  regression  rate  itself  is  continuous  at  the  breakpoints  observed  in 
nitramine  combustion,  although  the  slope  is  not.  If  the  reaction  site  (the  melt 
layer)  disappeared  completely,  we  would  expect  a discontinuity  in  the  regression 
rate,  as  indeed  the  AP  regression  rates  are  discontinuous  at  such  a point. 
However,  none  is  apparent  in  nitramine  combustion.  Similarly,  discontinuous 
variations  in  either  the  wall  temperature  or  the  fundamental  rate  limiting 
reaction  would  also  be  expected  to  be  reflected  in  discontinuities  in  the  regression 
rate,  except  in  the  fortuitous  circumstance  that  two  or  more  variations  occur  in 
a mutually  compensating  manner.  These  points  reinforce  the  argument  that  the 
basic  mechanisms  behind  the  regression  are  probably  not  too  different  trom 
each  other  on  e;ther  side  of  the  breakpoint. 

If  chemical  processes  were  of  crucial  importance  to  the  slope-break 
phenomenon,  it  would  seem  reasonable  to  expect  that  pressure  would  have  a 
much  stronger  effect  than  it  is  observed  to  have.  Kor  example,  we  ought  not 
to  be  able  to  shift  the  breakpoint  position  on  the  pressure  scale,  not  to  mention 
the  complete  elimination  of  the  breakpoint  (from  the  desired  operating  range,  at 


least)  through  variations  in  a physical  parameter  such  as  the  oxidizer  particle 


If  the  high  slope  is  caused  by  increased  surface  area  due  to  mechanical 
cracking  of  the  crystals,  one  would  expect  a fairly  random  behavior  of 
regression  rate  with  pressure  in  the  post-breakpoint  regime.  None  is  evident. 
The  data  seem  to  be  very  reproducible,  besides  showing  a well  behaved  n value. 

The  above  arguments  are  presented  to  indicate  the  strong  need  for  more 
experimental  data  and  are  not  intended  as  solutions  to  the  problem. 
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Section  v 


CONCLUDING  REMARKS 


A.  GENERAL  REMARKS 

The  present  theoretical  work  on  nitramine  propellant  combustion  has 
considered  a parametric  study  of  the  various  quantities  of  interest.  The 
ultimate  goal  of  all  such  studies  would  be  to  analytically  predict  the  regression 
rate  as  a function  of  the  propellant  formulation  and  the  chamber  pressure.  If 
successful  at  this  task,  we  would  be  in  a position  to  alter  the  formulations  to 
suit  specific  needs  of  a propellant  pressure-time  history.  In  fact,  that  was 
precisely  the  main  motivation  behind  the  present  investigations.  Since  we  have 
not  achieved  such  a predictive  ability  even  in  the  vastly  explored  AP/composites 
field,  it  is  needless  to  add  that  the  main  task  still  remains  open  in  the  nitramine 
propellant  field.  However,  this  work  has  indicated  the  general  trends  fairly 
well.  Probably  the  single  most  important  result  is  the  clear  indication  that  our 
model,  which  started  out  with  postulates  and  hypotheses,  is  indeed  applicable 
to  the  problem  of  nitramine  combustion.  The  fact  that  the  fundamental  rate  data 
used,  with  no  free  parameters,  predicts  the  observed  trends  cannot  be  indicative 
of  the  contrary. 

In  general,  our  model  has  relied  rather  heavily  on  AP/composite  pro- 
pellant data.  This  is  true  not  so  much  with  regard  to  the  actual  use  of  that 
data  but  more  with  regard  to  the  position  of  the  AP/propellant  field.  More 
specifically,  in  earlier  work,  the  AP/propellant  data  were  repeatedly  used 
(Refs.  2-4)  as  the  standard  testing  basis  for  our  model.  The  model  consis- 
tently predicted  results  close  to  experimental  observations.  Generalizations 
have  since  been  attempted  to  cover  the  nitramine  propellant  field.  It  is  natural 
to  ask  about  th  i validity  of  such  generalizations. 

The  similarities  are  thought  to  be  profound.  Both  of  these  propellants 
(AP  and  nitramine)  are  composite  propellants  using  a rather  heavy  loading  of 
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crystalline  oxidizers.  The  binders  form  a rather  small  portion  of  the  formulation 
and  usually  are  also  polymers  of  some  form  or  other.  The  degradation  reactions 
of  the  oxidizer  and  the  binder  seem  to  obey  the  same  general  type  of  Arrhenius 
expressions.  The  pre -exponential  factor  was  assumed  to  be  linearly  dependent 
on  the  chamber  pressure,  based  on  a picture  of  degradation.  Later  it  was 
found  that  experimenters  had  indeed  encountered  this  pressure-dependent 
degradation  of  AP  (Ref.  31).  Based  on  the  same  picture,  pressure-dependent 
degradation  of  nitramines  has  been  assumed  here.  We  are  hopeful  that  a future 
effort,  or  a future  discovery  of  an  earlier  experimental  study,  might  indicate 
the  assumed  pressure  dependence. 

An  extremely  important  aspect  of  the  present  study  needs  to  be  clearly 
understood.  Our  model  has  been  put  forward  to  represent  composite  propellants 
as  a family.  Except  for  the  numerical  constants,  the  model  makes  absolutely  no 
distinction,  at  least  in  its  present  form,  between  AP  composites  and  nitramine 
composites.  The  model  says  that  beyond  the  homogeneous  solid  limit  for 
composite  propellants  (i.e.,  when  a>*/r),  we  should  observe  changes  in  the 
propellant  combustion  behavior,  no  matter  what  the  chemical  lormulations  are. 
The  natural  question  would  be  to  inquire  whether  such  slope-break  or  similar 
phenomena  are  known  to  occur  with  AP/composite  propellants.  At  first  sight 
the  AP/propellant  family  may  appear  to  be  totally  devoid  of  all  such  special 
behavior.  However,  a closer  examination  of  the  data  indicates  that  most  AP/ 
propellants  have  been  meant  generally  for  rocket  applications  where  the 
operating  pressures  are  rather  low  compared  to  the  high  pressure  applications 
that  the  nitramine  propellant  data  cover.  Hence  we  should  examine  the  AP/ 
propellant  data  in  high  pressure  applications.  Fortunately,  such  a study  has 
been  made  and  the  report  has  been  declassified  recently  (Ref.  32).  That  study, 
which  involved  closed  bomb  tests  of  high  pressure  combustion  of  AP/composite 
propellants  of  a variety  of  formulations,  invariably  found  an  increase  in  the 
slcpe,  n,  by  almost  a factor  of  two  at  the  higher  regression  rate.  The  data 
presented  (Refs.  32-34)  is  indeed  very  impressive  and  shows  unmistakahle 
changes  in  the  slopes  at  high  pressures.  The  slope  breaks  appear  to  be  rather 
smooth  and  not  as  abrupt  as  those  reported  for  nitramines.  This  could  be  due 
to  the  multimodal  distribution  of  particles  in  the  conventional  AP/composites 
or  could  also  be  due  to  broader  distribution  of  particle  sizes  in  the  AP  case  as 


compared  with  the  nitramines  case.  [As  a matter  of  fact,  a fairly  sharp  break 
was  observed  (Ref.  26)  with  a fa.  rly  narrow  AP  size  distribution  in  a nonmetaLlized 
AP/PBAN  propellant.  ] Also,  when  an  extremely  unimodal  distribution  of  nitra- 
mine  oxidizer  crystals  was  used  in  a propellant  (Ref.  15),  a drastically  high 
slope  break  was  discovered.  The  main  ’’signal"  in  all  of  this  is  simply  that  the 
slope  change  anticipated  by  the  CIT/JPL  model  for  composite  propellants  in 
general  is  indeed  observed  not  only  in  nitramine  propellants,  but  in  the  conven- 
tional AP/composites  as  well.  This  must  certainly  be  regarded  as  a strong 
support  for  our  model. 

The  second  aspect  of  this  work  has  to  do  with  specific  propellant 
applications.  All  of  the  numerical  data  were  obtained  for  a model  composite 
propellant  that  has  the  properties  of  RDX  oxidizer.  Rate  constants  were  found 
(Ref.  12)  only  for  RDX.  In  any  case  the  results  are  not  expected  to  be  widely 
different  for  the  other  nitramine  oxidizers.  Superficially,  it  may  appear  that 
binder  interactions  have  been  completely  ignored  in  this  work.  Actually,  it  is 
the  binder  interaction  that  determines  such  important  parameters  as  the  wall 
temperature  T , the  flame  temperature  T^,  and  the  other  property  values  such 
as  p,  c,  k , etc.  Thus,  the  binder  interactions  have  been  intimately  interwoven 
into  our  model.  Also,  the  gas-phase  model  used  for  AP/composites  may  appear, 
at  first  sight,  to  be  totally  inapplicable  to  the  nitramine  composites.  This  is 
mainly  because  of  the  fundamentally  diflerent  roles  that  the  binder  vapors  play 
in  the  combustion  of  these  two  classes  of  propellants.  In  AP/composites  the 
AP  monopropellant  flame  temperature  is  lower  than  the  propellant  flame 
temperature;  the  binder  vapors  undergo  chemical  reactions  with  the  oxidizer 
vapors  exothernially  to  reach  the  final  flame  temperature.  In  the  nitramine 
propellants  with  the  so-called  "inert"  binders,  the  monopropellant  oxidizer 
tlame  temperature  is  higher  than  the  propellant  flame  temperature  and  the 
binder  vapors  actually  lower  the  temperature  with  interaction.  However,  under 
the  basic  assumption  that  the  diffusion-mixing  process  of  the  two  vapors 
(oxidizer's  and  binder's)  controls  the  gas-ohase  temperature  profile,  the  same 
vapor-phase  model  is  applicable  for  both  the  propellants.  In  addition,  we  have 
the  interesting  situation  that  the  same  numerical  constants  used  for  AP/compo- 
sites  may  be  used,  as  a first  approximation,  in  describing  this  mixing  process 
(and  hence  the  heat  release  process).  The  important  mixing  process  is 
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nonlaminar  and  this  physical  diffusion-mixing  is  unlikeLy  to  depend  strongly  on 
the  chemical  nature  of  the  vapors.  In  any  case,  the  gas-phase  processes  are 
not  the  primary  controlling  factors  in  propellant  burning  in  our  model. 

It  is  encouraging  that  the  predictions  and  the  anticipations  of  one  model 
are  consistent  with  the  main  trends  in  the  experimental  literature.  The  model 
indicates  that  decreasing  the  oxidizer  particle  size  should  have  a beneficial 
effect;  it  is  indeed  observed  to  have  a beneficial  effect.  The  model  places 
little  emphasis  on  the  metal  content,  special  ingredients,  and  the  chemical 
nature  of  the  oxidizer.  Experimentally,  it  is  found  that  "whether  inert  or 
energetic  binder,  cyclic  or  linear  nitramine,  metallized  or  non-metallized,  the 
slope-break  phenomenon  is  invariably  exhibited"  (Ref.  1). 

In  the  context  of  our  model,  the  catalyst  would  be  influencing  tin  ondensed- 
phase  degradation  rate  and  consequently  the  break  in  the  regression  rate  vs. 
pressure  would  persist  around  the  same  value  of  tne  regress. on  rate.  However, 
the  pressure  at  the  breakpoint  would  be  lower,  because  the  catalyst  addition 
would  increase  the  regression  rate  at  a given  pressure.  Hence,  in  the  context 
of  the  CIT/JPL  model,  it  would  appear  that  the  breakpoint  would  occur  at  lower 
pressures  with  catalyst  addition,  if  no  other  changes  are  made  at  the  same  time. 
{However,  the  catalyst  itself  is  obviously  added  at  the  expense  of  some  other 
component,  and  we  do  not  have,  strictly  speaking,  the  same  composition  any 
more. ) 

Our  model,  admittedly,  is  highly  idealized  and  needs  improvements  before 
handling  more  realistic  features  such  as  particle  size  distributions,  multimodal 
oxidizers,  etc.  Also,  further  studies  are  expected  to  make  the  model  more  self- 
contained  in  that  the  wall  temperature  variations  will  be  automatically  handled. 

To  summarize  those  aspects  that  affect  the  design  of  propellants,  our  model 
indicates  that  the  following  steps  ought  to  help  obviate  the  high  slope  from  the 
operating  range: 
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1)  Decreasing  the  oxidizer  particle  size. 

2)  Employing  a binder  system  that  melts  readily  (PU,  for  example). 

4)  Increasing  the  thermal  diffusivity  of  the  propellant. 

4)  In  general,  introducing  any  modification  that  enhances  the 

homogeneity  of  the  condensed  phase. 

It  is  believed  that  mere  research  effort  will  actually  enable  us  to  arrive 
theoretically  at  propellant  formulations  to  meet  specific  applications  criteria. 

B.  SALIENT  CONCLUSIONS 

The  above  parametric  study  on  the  combustion  of  nitramine  propellants 
has  led  to  the  following  conclusions: 

1)  The  CIT/JPL  model  realistically  predicts  the  observed  trends 
with  a minimum  number  of  input  parameters. 

2)  Different  gas-phase  models  lead  us  to  essentially  the  same 
analytical  trends,  thus  de -emphasizing  the  importance  of  the 
gas-phase  details. 

3)  The  wall  temperature  of  the  propellant  decreases  with  increases 
in  oxidizer  particle  size;  the  interfacial  melt  layer  thickness  has 

a relatively  minor  influence  on  the  propellant  combustion  behavior. 

4)  The  hypothesis  of  gas-phase  combustion  zone  dimensions  being 
determined  by  the  condensed-phase  heterogeneity,  a,  leads  to 
predictions  of  the  wall  temperature  ana  the  regression  rate 
that  are  both  in  agreement  with  the  experimentally  observed 
trends . 

5)  With  regard  to  the  highly  popular  belief  in  the  existence  of  a "melt 
layer"  on  the  surface  of  the  propellant,  the  model  seems  to  be 
sufficiently  general  to  be  consistent  with  the  reported  observations. 

f>)  Simplified  relationships  for  the  melt  layer  thickness  and  flame 

standoff  distance  give  reasonable  results  which  show  them  decreasing 
with  increasing  pressures,  consistent  with  reported  observations. 
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This  report  summarizes  a theoretical  investigation  of  the  current  problems  of 
nitmmine  (composite)  propellant  combustion.  This  study  has,  as  its  distinctive 
featuie,  a detailed  examination  of  the  condensed-phase  processes  in  the  combustion 
cf  nitramine  pi  ope  it  ants.  As  a consequence  of  a recently  developed  model  for  the 
combustion  of  anunonMin  perchlorate  (AP)/composite  propellants,  it  is  hypothesized 
that  the  condensed-phase  degradation  of  the  nitramine  oxidizer  particles  to  a 
vaporizable  state  is  the  overall  rate-limiting  step.  It  is  also  assumed  that  the 
gas -phase  details  are  secondary  in  importance  and  need  be  studied  only  to  the 
extent  of  supplying  the  correct  boundary  conditions  on  the  condensed-phase/vapor- 
phase  heat  transfer.  Because  of  our  imprecise  understanding  of  the  gas-phase 
processed  in  the  presence  of  combustion,  several  plausible  models  are  considered 
or  the  gas  phase.  It  is  found  that  all  of  the  gas-phase  models  considered  lead 
to  predictions  sufficiently  close  to  experimental  trends  for  us  to  conclude  that 
the  precise  details  of  gas-phase  processes  are  not  of  critical  importance  in 
determining  propellant  combustion  behavior.  More  to  the  point,  we  are  led  to 
believe  that  a thorough  examination  of  the  condensed-phase  details  may  be  suffi- 
cient in  itself  not  only  to  interpret  most  of  the  available  data  on  experimental 
regression  rate  vs.  pressure  of  nitramine  propellants,  but  also  to  aid  in  the 
formulation  of  propellants  to  suit  our  needs. 
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